Regimes of quantum degeneracy in trapped ID gases 
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We discuss the regimes of quantum degeneracy in a trapped ID gas and obtain the diagram of states. 
Three regimes have been identified: the BEC regimes of a true condensate and quasicondensate, 
and the regime of a trapped gas of Tonks (gas of impenetrable bosons). The presence of a sharp 
cross-over to the BEC regime requires extremely small interaction between particles. We discuss 
how to distinguish between true and quasicondensates in phase coherence experiments. 
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Low-temperature ID Bose systems attract a great deal 
of interest as they show a remarkable physics not encoun- 
tered in 2D and 3D. In particular, the ID Bose gas with 
repulsive interparticle interaction (the coupling constant 
g > 0) becomes more non-ideal with decreasing ID den- 
sity n Q . The regime of a weakly interacting gas requires 
the correlation length l c = h/^/mng (m is the atom mass) 
to be much larger than the mean interparticle separation 
1/n. For small n or large interaction, where this condi- 
tion is violated, the gas acquires Fermi properties as the 
wavefunction strongly decreases at short interparticle dis- 
tances . In this case it is called a gas of impenetrable 
bosons or gas of Tonks (cf. 

Spatially homogeneous ID Bose gases with repulsive 
interparticle interaction have been extensively studied in 
the last decades. For the delta-functional interaction, 
Lieb and Liniger Q| have calculated the ground state en- 
ergy and the spectrum of elementary excitations which 
at low momenta turns out to be phonon-like. Generaliz- 
ing the Lieb-Liniger approach, Yang and Yang have 
proved the analyticity of thermodynamical functions at 
any finite temperature, which indicates the absence of 
a phase transition. However, at sufficiently low T the 
correlation properties of a ID Bose gas are qualitatively 
different from classical high-T properties. In the regime 
of a weakly interacting gas (nl c 3> 1) the density fluctua- 
tions are suppressed |5| , whereas at finite T the long- wave 
fluctuations of the phase lead to exponential decay of the 
one-particle density matrix at large distances HH- A 
similar picture, with a power-law decay of the density 
matrix, was found at T = [Rj8|. Therefore, the Bose- 
Einstein condensate is absent at any T, including T = 0. 
Earlier studies of ID Bose systems are reviewed in 
They allow us to conclude that in ID gases the decrease 
of temperature leads to a continuous transformation of 
correlation properties from ideal-gas classical to interac- 
tion/statistics dominated. A ID classical field model for 
calculating correlation functions in the conditions, where 
both the density and phase fluctuations are important, 
was developed in Q and with respect to Bose gases in 
jnj. Interestingly, ID gases can posses the property of 



superfluidity at T = jl2]j|]. Moreover, at finite T one 
can have metastable supercurrent states which decay on 
a time scale independent of the size of the system jll| . 

The earlier discussion of ID Bose gases was mostly 
academic as there was no possible realization of such a 
system. Fast progress in evaporative and optical cooling 
of trapped atoms and the observation of Bose-Einstein 
condensation (BEC) in trapped clouds of alkali atoms 
|l4| stimulated a search for non-trivial trapping geome- 
tries. At present, there are significant efforts to create 
(quasi) ID trapped gases where the radial motion of 
atoms in a cylindrical trap is tightly confined and they 
only undergo zero point radial oscillations. Then, kine- 
matically the gas is ID, and the difference from purely 
ID gases will only be related to the value of the interpar- 
ticle interaction which now depends on the radial con- 
finement. The presence of the axial confinement allows 
one to speak of a trapped ID gas. 

Ketterle and van Druten Jl6| considered a trapped ID 
ideal gas and have revealed an essential role of the dis- 
crete structure of the trap levels. They have established 
that at temperatures T < Nhui / In 2N, where N is the 
number of particles and lo the trap frequency, the popu- 
lation of the ground state rapidly grows with decreasing 
T and becomes macroscopic. Thus, one has a clear BEC- 
like behavior of the ideal trapped ID cloud. 

A fundamental question concerns the influence of in- 
terparticle interaction on the presence and nature of a 
Bose-condensed state and on the character of a cross-over 
to the BEC regime. In this Letter we discuss the regimes 
of quantum degeneracy in a trapped ID gas with repul- 
sive interparticle interaction. We find that the presence 
of a sharp cross-over to the BEC regime, predicted in 
JLsf , requires extremely small interaction between parti- 
cles. Otherwise, the decrease of temperature leads to a 
continuous transformation of a classical gas to quantum 
degenerate. We identify 3 regimes at T < Td, where 
Td rs Nfobj is the degeneracy temperature. For a suffi- 
ciently large interparticle interaction and the number of 
particles much smaller than a characteristic value , at 
any T -C Td one has a trapped gas of Tonks, with the 
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density profile characteristic for an ideal Fermi gas. For 
N ^> N* we have a weakly interacting gas. The pres- 
ence of the trapping potential introduces a finite size of 
the sample and drastically changes the picture of long- 
wave fluctuations of the phase compared to the earlier 
discussed uniform case. We calculate the density and 
phase fluctuations and find that well below Td there is 
a quasicondensate, i.e. the Bose-condensed state where 
the density fluctuations are suppressed but the phase still 
fluctuates. At very low T the long-wave fluctuations of 
the phase are suppressed due to a finite size of the system, 
and we have a true condensate. The true condensate and 
the quasicondensate have the same Thomas-Fermi den- 
sity profile and local correlation properties, and we an- 
alyze how to distinguish between these Bose-condensed 
states in phase coherence experiments. 

We first discuss the coupling constant g for possible 
realizations of ID gases. These realizations imply parti- 
cles in a cylindrical trap, which are tightly confined in 
the radial (p) direction, with the confinement frequency 
ujq greatly exceeding the mean-field interaction. Then, at 
sufficiently low T the radial motion of particles is essen- 
tially "frozen" and is governed by the ground-state wave- 
function of the radial harmonic oscillator. If the radial 
extension of the wavefunction, Iq — (h/mwo) 1 ^ 2 ^> R e , 
where R e is the characteristic radius of the interatomic 
potential, the interaction between particles acquires a 3D 
character and will be characterized by the 3D scattering 
length a. In this case, assuming Iq 3> \a\, we have 



g = 2h 2 a/mlQ. 



(1) 



This result follows from the analysis in |17j and can also 
be obtained by averaging the 3D interaction over the ra- 
dial density profile. Thus, statistical properties of the 
sample are the same as those of a purely ID system with 
the coupling constant given by Eq.(|l]). 

In the regime of a weakly interacting gas, where nl c 3> 
1, we have a small parameter 



7 = l/(nl c ) 2 = mg/h 2 n < 1. 



(2) 



For particles trapped in a harmonic (axial) potential 
V(z) = muj 2 z 2 /2, one can introduce a complementary 
dimensionless quantity a — mgl/h 2 which provides a re- 
lation between the interaction strength g and the trap 
frequency ui (I = y'R/mw is the amplitude of axial zero 
point oscillations). 

At T = one has a true condensate: In the Thomas- 
Fermi (TF) regime the mean square fluctuations of the 
phase do not exceed ~ 7 1 / 2 and, hence, these fluctuations 
are small under the condition (||) (see [18| ) . The conden- 
sate wavefunction is determined by the Gross-Pitaevskii 
equation (GPE) which gives the well-known TF parabolic 
density profile no(z) = riQ ln (l — z 2 /R 2 ^ F ). The maximum 
density no m = n(0) = fi/g, the TF size of the conden- 
sate Rtf — (2fi/muj 2 ) 1 ^ 2 , and the chemical potential 



(i — huj(3Na/4V2) 2 / 3 . For a > 1 we are always in the 
TF regime (p huj). In this case, Eq.(||) requires a 
sufficiently large number of particles: 



N 3> = a 



(3) 



Note that under this condition the ratio fi/Td ~ 
(a 2 /N) 1/3 < 1. For a < 1 the criterion (j|) of a weakly 
interacting gas is satisfied at any N, and the condensate 
is in the TF regime if N ^> a -1 . In the opposite limit 
the mean-field interaction is much smaller than the level 
spacing in the trap hu>. Hence, one has a macroscopic 
occupation of the ground state of the trap, i.e. there is 
an ideal gas condensate with a Gaussian density profile. 

At this point, we briefly discuss the cross-over to the 
BEC regime, predicted by Ketterle and van Druten [fl6|| . 
They considered a trapped ID ideal gas and found that 
the decrease of temperature to below T c = NTilo / In 2N 
leads to a strong increase of the population of the ground 
state, which rapidly becomes macroscopic. This sharp 
cross-over originates from the discrete structure of the 
trap levels and is not observed in quasiclassical calcula- 
tions |19|. We argue that the presence of the interparticle 
interaction changes the picture drastically. One can only 
distinguish between the (lowest) trap levels if the inter- 
action between particles occupying a particular level is 
much smaller than the level spacing. Otherwise the in- 
terparticle interaction smears out the discrete structure 
of the levels. For T close to T c the occupation of the 
ground state is ~ T c /Tilo « N/h\2N |ll| and, hence, the 
mean-field interaction between the particles in this state 
(per particle) will be Ng/l\n2N. The sharp BEC cross- 
over requires this quantity to be much smaller than fiu), 
and we arrive at the condition N/ In 2 N <C a -1 . For a 
realistic number of trapped particles (N ~ 10 3 — 1 4 ) 
this is practically equivalent to the condition at which 
one has the ideal gas Gaussian condensate (see above). 

As we see, the sharp BEC cross-over requires small a. 
For possible realizations of ID gases, using the coupling 
constant g (|l|), we obtain a — 2o,1/Iq. Then, even for the 
ratio I/Iq ~ 10 and moderate radial confinement with 
Iq ~ 1/zm, we have a ~ 0.1 for Rb atoms (a w 50 A). 
Clearly, for a reasonably large number of particles the 
cross-over condition JV < a" 1 can only be fulfilled at 
extremely small interparticle interaction. One can think 
of reducing a to below lA and achieving a < 10 -3 by 
using Feshbach resonances as in the MIT and JILA ex- 
periments (2^]. In this case, already for N ~ 10 3 one can 
expect the sharp BEC cross-over and the existence of the 
Gaussian condensate at T <T C . 

We now turn to the case of large a. For 7 ^> 1 
one has a gas of Tonks ]ipl|]. The one-to-one map- 
ping of this system to a gas of free fermions || en- 
sures the fermionic spectrum and density profile of a 
trapped gas of Tonks. For (axial) harmonic trapping 
the condition 7 3> 1 requires iV < JV,. The chemi- 
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cal potential is equal to Nfruj, and the density distribu- 
tion is n(z) = n m yjl — (z/R) 2 , where n m = ^/2N/nl, 

f2NL The density pro- 



n m ^l - (z/R) 2 , 
and the size of the cloud R = \ 
file n(z) is quite different from both the profile of the 
zero-temperature condensate and the classical distribu- 
tion of particles, which provides a root for identifying 
the trapped gas of Tonks in future experiments. The in- 
terference effects and dynamical properties of this system 
are now a subject of theoretical studies ]22|j23| ]. In Rb and 
Na this regime can be achieved for TV < 10 3 by the Fes- 
hbach increase of a to ~ 500 A and using u> ~ 1Hz and 
optical radial confinement with ujq ~ 10kHz (a ~ 50). 

Large a and N satisfying Eq.(0), or small a and N » 
a , seem most feasible in experiments with trapped ID 
gases. In this case, at any T <C one has a weakly in- 
teracting gas in the TF regime. Similarly to the uniform 
ID case, the decrease of temperature to below T<j con- 
tinuously transforms a classical ID gas to the regime of 
quantum degeneracy. At T = this weakly interacting 
gas turns to the true TF condensate (see above). It is 
then subtle to understand how the correlation properties 
change with temperature at T <C Tj. For this purpose, 
we analyze the behavior of the one-particle density ma- 
trix by calculating the fluctuations of the density and 
phase. We a priori assume small density fluctuations 
and prove this statement relying on the zero-temperature 
equations for the mean density Uq(z) and excitations. 
The operator of the density fluctuations is (see, e.g. 
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where a,j is the annihilation operator of the excitation 
with quantum number j and energy ej, /• = Uj ± Vj, 
and the u, v functions of the excitations arc determined 
by the same Bogolyubov-de Gennes equations as in the 
presence of the TF condensate. 

The solution of these equations gives the spectrum 6j — 



hujy/ j(j + l)/2 [18 24 1 and the wavefunctions 
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where j is a positive integer, Pj are Legendre polynomi- 
als, and x = z/ Rtf- For the mean square fluctuations of 
the density, ((Sn(z, z')) 2 ) = ((h'(z) — h'(z')) 2 ), we have 



{(5n(z,z')f 



L 0m 



E 



; -(j + l/2) 
^ 2fJ.n 0m RTF 



(P j (x)-P j (x')) 2 (l + 2N j ) 



with Nj = [exp(ej/T) — being the occupation num- 
bers for the excitations. At T > hu the main contribu- 
tion to the density fluctuations comes from quasiclassical 
excitations (j ^> 1). The vacuum fluctuations are small: 
((Sn(z, z')) 2 )q ~ tiq^ 1 / 2 . For the thermal fluctuations 
on a distance scale \z — z'\ 3> l c , we obtain 



{(5n(z,z')f 
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We see that the density fluctuations are strongly sup- 
pressed at temperatures T <C Then, one can write the 
total field operator as ip(z) — y/no(z) exp(i(f>(z)) , where 
4>(z) is the operator of the phase fluctuations, and the 
one-particle density matrix takes the form (see, e.g. Pj) 



(^(z)^( Z ')}^y/n (z)n (z')eM-(mz,z')} 2 )/2}, (7) 

with 5<fi(z, z') — <fi(z) — 4>(z'). The operator of the phase 
fluctuations is given by (see p5|) 



4>(z) = [An Q (z)]- 1 / 2 Y,f^(z)a J +h.c, 

3=1 

and for the mean square fluctuations we have 



(8) 
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For the vacuum fluctuations we find (c.f. j26| . [18[ ) 

([^(z,z')] 2 )o«(7 1/2 A)ln(k-^'|/g, 

and they are small for any realistic size of the gas cloud. 
The thermal fluctuations of the phase are mostly pro- 
vided by the contribution of the lowest excitations. A 
direct calculation, with Nj — T/ej, yields 



([^(z,z')] 2 )t = 



4T/i 
3Tdhto 



log 



(1 - x') (1 + x) 



(1 



(1-x) 



(9) 



In the inner part of the gas sample the logarithm in 
Eq. (j^) is of order unity. 

Thus, we can introduce a characteristic temperature 



T p h = T d hw/fj, 



(10) 



at which the quantity ([64>(z, z')] 2 ) s» 1 on a distance 
scale \z — z'\ ~ Rtf- The characteristic radius of phase 
fluctuations is R^ « R TF (T ph /T) cx iV 2 / 3 /T, and for 
T < Tph it exceeds the sample size Rt f ■ This means that 
at T <C Tph both the density and phase fluctuations are 
suppressed, and there is a true condensate. The condition 
@ always provides the ratio T ph /huj w (AN /a 2 ) 1 / 3 > 1. 

In the temperature range, where Td ^> T 3> T p h, 
the density fluctuations are suppressed, but the phase 
fluctuates on a distance scale R^ -C Rtf- Thus, simi- 
larly to the quasi2D case p7|] , we have a condensate with 
fluctuating phase (quasicondensate) . The radius of the 
phase fluctuations greatly exceeds the correlation length: 
R$ ~ l c (Td/T) 3> l c . Hence, the quasicondensate has 
the same density profile as the true condensate. Correla- 
tion properties at distances smaller than R^ are also the 
same. However, the phase fluctuations lead to a drastic 
difference in the phase coherence properties. 
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In Fig.l. we present the state diagram of the trapped 
ID gas for a = 10 (A* = 100). For A > A*, the decrease 
of temperature to below Tj, leads to the appearance of 
a quasicondensate which at T < T p h turns to the true 
condensate. In the T — A plane the approximate border 
line between the two BEC regimes is determined by the 
equation {T/Tiuj) = (32A/9A,) 1 / 3 . For N < A» the 
system can be regarded as a trapped gas of Tonks. 



a trapped gas of Tonks. The creation of ID gases will 
open handles on interesting phase coherence studies and 
the studies of "fcrmionization" in Bose systems. 
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FIG. 1. Diagram of states for a trapped ID gas. 

Phase coherence properties of trapped ID gases can be 
studied in 'juggling' experiments similar to those with 3D 
condensates at NIST and Munich Small clouds 

of atoms are ejected from the main cloud by stimulated 
Raman or RF transitions. Observing the interference 
between two clouds, simultaneously ejected from differ- 
ent parts of the sample, allows the reconstruction of the 
spatial phase correlation properties. Similarly, temporal 
correlations of the phase can be studied by overlapping 
clouds ejected at different times from the same part of 
the sample. In this way juggling experiments provide a 
direct measurement of the one-particle density matrix. 
Repeatedly juggling clouds of a small volume O from 
points z and z' of the sample, for equal time of flight 
to the detector we have the averaged detection signal 
/ = fi[noC0 + n (z') + 2(ft(z)TP(z'))\. 

At T <C Tph the phase fluctuations are small and one 
has a true condensate. In this case, for z' = —z we 
have (z)4>(z')) = n (z) and / — 4nn (z), and there 
is a pronounced interference effect: The detected signal 
is twice as large as the number of atoms in the ejected 
clouds. The phase fluctuations grow with T and for 
T > Tph, where the true condensate turns to a quasi- 
condensate, the detection signal decreases as described 
by (^(z)^(z')) from Eqs. (@) and (§. For T > T ph 
the phase fluctuations completely destroy the interfer- 
ence between the two ejected clouds, and I — 2ilno(z). 

In conclusion, we have identified 3 regimes of quantum 
degeneracy in a trapped ID gas: the BEC regimes of a 
quasicondensate and true condensate, and the regime of 
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